High-density correlation energy expansion of the one-dimensional uniform electron gas 
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We show that the expression of the high-density (i.e small-r 3 ) correlation energy per electron for 
the one-dimensional uniform electron gas can be obtained by conventional perturbation theory and 
is of the form e c (r s ) = — 7r 2 /360 + 0.00714r s + . . ., where r s is the average radius of an electron. 
Combining these new results with the low-density correlation energy expansion and the available 
diffusion Monte Carlo (DMC) data, we propose a local-density approximation correlation functional, 
which deviates by a maximum of 0.2 millihartree compared to the benchmark DMC calculations. 
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I. INTRODUCTION 

While the high-density (small-r s ) reduced (i.e. per electron) correlation energy expansions 

oo 

e c (V s ) = ( A j lnr » + e j) r l 

j=o I 1 ) 

= A In r s + e + Xir s In r s + eir s + . . . 

(where r s is the Seitz radius) of the two-dimensional (2D) and three-dimensional (3D) uniform electron gas (UEG) are 
quite well-knownpN23 much less has been discovered about the one-dimensional (ID) UEG. This lack of information 
is mainly due to the divergence of the Coulomb operator 1/x in ID for small interelectronic distance xp^^ which 
makes conventional perturbation theory difficult to apply due to the absence of Fourier transform for the Coulomb 
operator. Recentl y m uch attention has been devoted to such ID systems, which can be experimentally realized in 
carbon nanotubesp^H organic conductors J2£H2ID transition metal oxides,Bll edge states in quantum Hall li quids BUBS 
semiconductor heterostructuresBSBH confined atomic gasesJ^HSU anc | atomic or semiconducting nanowires ! 53 * 54 * 

In this Article, we propose to fill this gap by reporting the values of the first few high-density coefficients (see Table 
HI . The present system is constructed by allowing the number n of electrons in a ID box of length L with periodic 
boundary conditions to approach infinity with the density p — n/L held constant ) 55 l 56 l Because the paramagnetic and 
ferromagnetic states are degenerate for strict ID systems, we will consider only the latterPStH] 

To avoid the divergence of the Coulomb operator, we will consider in our derivation a "sof ten" version of the 
Coulomb operator l/\/x 2 + R 2 , where R is a parameter which removes the singularity at x = o) 57 l 58 l Then, we will 
carefully take the limit R — > 0. We will show that, unlike the 2D and 3D version of the UEG, second- and third-order 
perturbation theories are convergent, i.e there is no need to use resummation techniques.^ Combining these new 
results with the low-density energy expansion and the available diffusion Monte Carlo (DMC) data, we propose a new 
local-density approximation (LDA) functional for the reduced correlation energy of the ID UEG. Atomic units are 
used throughout. 



II. HIGH-DENSITY EXPANSION 
A. Second-order perturbation theory 



In ID, the spinorbitals of the free electron gas are 



M*) = -7^' ( 2 ) 



with the energy = k 2 /2, and where the periodic boundary conditions imply k = 2irm/L (m 6 Z). The coefficient 
eo is given by second-order perturbation theory^ 



art ' * ' * k - -A- K.%. — /c„ — k. _ 

where (a6||rs) = (ab\rs) — (ab\sr) and 



ab rs 



{ a b \rs) = r r ^;^^>^Wdx 2 . (4) 

J-L/2J-L/2 - x 2 ) 2 + R 2 

The constant coefficient eo is usually decomposed in a direct ("ring-diagram") term €q and an exchange term Cq, which 
read explicitly 



occ virt 
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TABLE I. High-density coefficients for the paramagnetic state of the ID, 2D and 3D UEGs. The paramagnetic and ferromagnetic 
states are degenerate in ID. j3 and £ are the Dirichlet beta and Riemann zeta function, respectively™ 



Coefficient 



Term 



ID 



2D 



3D 



Ac 
Ai 

61 



In r 3 

r° 3 -tt 2 /360 In 2 - 1 + /9(2) - 8/3(4) /V 2 

r s \nr s -y/2 (10/(3ir) - 1) 

r s +0.00714 unknown 



(1 -In2)/7T 2 
-0.071099 + (ln2)/6 - 3C(3)/(4tt 2 ) 
+0.009229 
-0.020 
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FIG. 1. Convergence of eo and ei with respect to R. 



Using the Fourier transform of the soft Coulomb potential 

1 1 



K {\k\Ry kx dk, (7) 



VX 2 + R 2 7T 

where Kq is the zeroth-order modified Bcssel function of the second kind™, the well-known relation 



6(x) = -j= / e ik *dk, (8) 

V27T J-oo 

and transforming the sums in ^ and (|6j) into integrals with p\ = a/k-p, P2 — b/kp and q = fe/fcp, where fcp = 717? is 
the Fermi wave vector, we eventually find^ 

° V ; ^ 2 J-oc J t ^<\ J t ^<\ q(Pi+P2 + q) 



and 

|Pl+<j|>l |P2+<?|>1 

For R > 0, Eqs. ((9j) and ((lOj can be evaluated numerically. As shown in Fig. [I] eo(R) decreases monotonically to 
reach a constant at R = 0. However, for R = 0, both integral diverge with opposite rate. Thus, to find the limiting 
value, it is better not to split eo{R) into two contributions but to consider them together. For small i?, we have 

|pi+?l>l b2+g|>i 
The integrations over p\ and P2 can be performed at this stage and it yields 

€ °( R ) = ~i ^pH<Z^/2) + 7]*Z, (12) 
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where 7 is the Euler-Mascheroni constant)^ and 

A(g) = (2 + g)[ln(2 + q)] 2 + (2 - g)[ln(2 - q)] 2 
-2(l + lng)(2 + g) ln(2 + g) 
-2(l + lng)(2-g) ln(2-g) 
-4(ln2) 2 + 81n2(l + ln g ) 

for < g < 2, and 

A(g) = (2 + g)[ln(2 + q)] 2 + (2 - g)[ln(g - 2)] 2 
-2(l + lng)(2 + g) ln(2 + g) 
-2(l + lng)(2-g) ln(g - 2) 
+ 2(2 + \nq)q\nq 



(13) 



(14) 



otherwise. Performing the last integration over the two distinct regions (0 < q < 2 and q > 2) gives two contributions 
which diverges as lni? for small R with opposite sign. Thus, the divergences cancel and we find 



n 2 



eo=limeo(E) = --, (15) 

which nicely reproduces the result obtained with a ring geometry.^ 

Because e c (r s ) = eo + 0(r s ) (see below), eo provides the exact value of the correlation energy at r s = 0, and it 
is roughly —27.4 millihartree per electron. It is worth noting that, in most of the studies on ID systems, a soft 
Coulomb operator is considered. For example, in Refs. W2\ and 1531 the authors used R = 1, yielding a correlation 
energy (-8.7 millihartree) more than three times lower than the value obtained using the genuine Coulomb operator (i.e 
R = 0). Moreover, using a quasi-lD model with a transverse harmonic potential, Casula et al. conclude that, in the 
high-density limit, the correlation energy vanishes quadratically with r s !^This strikingly different prediction stresses 
the importance of employing a realistic Coulomb operator 



B. Third-order perturbation theory 

Using the same approach, third-order perturbation theory giveiP^ 

ei( 



1 ^ (ab\\rs)(cd\\ab)(rs\\cd) 



8n \K a + Kf) K r K s )(k c + Krf K r K s ) 

abcdrs 



1 v (ab\\rs)(rs\\tu}(tu\\ab) 

8n f-^ (n a + n b - K r - K s )(n a + n b - K t ~ k u ) ( 16 <* 

abrstu 

1 ^ (ab\\rs)(cs\\tb)(rt\\ac) 



^ abcrst ^ a ^ ^ r ^s){^a ~\~ ^t) 



Equation (16) can be decomposed, using the same transformations as in ^ and into six distinct contributions: 

7r 4 J q{pi 



+ P2 + q)q(pi +P3 + q) 

, ,3 v R (qfv R (pi -p 2 ) 
dq d°p— — r 

q{p\ + ps + q)q{P2 + P3 + q) 

2 f dqd 3 p- Vr ^ 2vr ^ 1 + p2 + q "> 



7T 4 J q(pi +P2 + q)q{pi + P3 + q) ^ 
2 f d d3 v R {q)v R (pi +p 2 + q)v R {p 2 - Pa) 



7T 4 J q(pi + P2 + q)q(pi + P3 + q) 

J_ f d 2 d 2 VR(qi) v n(q2) [vR(qi - ga) - vr(pi - P2)] 
7T 4 J qi(pi -P2)?2(pi -P2) 

J_ f d 2 d 2 v R(qi)vR(q2) [vR(qi - gg) - vr(pi +P2 + gi + g 2 )] 

27r 4 J q 1 (p 1 +p 2 + <7i)<72(pi +P2 + 92) 
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TABLE II. Reduced correlation energy (—e c (r 3 ) in millihartree) for various r s . The DMC and ISI results are taken from 
Refs. l28l and 1301 respectively. Subscripts represent the statistical errors in the last digits. The deviation with respect to the 
DMC result is given in parenthesis. 





DMC 


ISI 


This work 


o 




27.416 


27.416 


0.1 




26.632 


26.714 


0.2 


— 


25.900 


26.039 


0.5 




23.964 


24.191 


1.0 


21.4441 2 


21.392 (-0.052) 


21.648 (+0.204) 


2.0 


17.922 02 7 


17.763 (-0.159) 


17.990 (+0.068) 


5.0 


12.31774a 


12.091 (-0.227) 


12.226 (-0.092) 


10.0 


8.292 096 9 


8.121 (-0.171) 


8.196 (-0.096) 


15.0 


6.319 404 4 


6.195 (-0.124) 


6.246 (-0.073) 


20.0 


5.132 504 2 


5.039 (-0.094) 


5.077 (-0.056) 




FIG. 2. e c (r s ) given by Eq. (21 1 as a function of r s (solid li ne). DMC results are shown by black dots. The small-r s expansion 



of Eq. ( 19 \ (dashed line) and large-r s approximation of Eq. ( 20 1 (dotted line) are also shown 



where vn{q) = K (\q\R) and / d m p = J ■ ■ ■ J dpi . . . dp m . The regions of integration are \pi\ < 1 and \pt + q\ > 1 for 
the first four integrals, and the fifth integral is over the region: \pi + qi\ > 1, \pi + cfel > 1, \p\\ < 1; \p2 + Qi\ < 1, 
\P2 + 92 1 < 1, |p2 1 > 1- The final integral is the sum of two parts, one over the region \pi + qi\ > 1, \pi + ^1 > 1, 
\pi \ < 1; \p2 + qi\ > 1, \p2 + 92 1 > 1, \P2\ < 1; the second over the region the opposite inequalities. 

Again, for R = 0, most of the integrals diverge. The first five terms have to be considered together, while the two 
parts of the last integral are finite. Evaluating numerically each of the contributions and extrapolating the result to 
R = using the relation aR@ + t\ (see Fig. nl), we find 

ei = Urn ei(.R) = 0.00714(1), (18) 



where the error has been obtained by taking into account each numerical error estimate and extrapolating the overall 
error to R ~ 0P^ W e note that the present ID UEG is one of the few systems where the r s coefficient of the high-density 
expansion is known ESESl Moreover, the "line" geometry considered here seems to be more mathematically convenient 
(at least numerically) for deriving the r s coefficient than the ring geometry of Ref. |30l 

In summary, we have shown that the high-density correlation energy expansion ([T]) of the ID UEG is 

e c (r s ) = -^r+ 0.00714 r 3 + . . . . (19) 



We note that, contrary to the 2D and 3D UEGs, the expansion ( 19 ) does not contain any logarithm term up to first 
order in r s , i.e. Aq — X\ = (cf Eq. ([T])). 
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III. LDA FUNCTIONAL 

It is knowrPSHH that the low-density (large-r s ) expansion of the correlation energy is 



"s r s 

(20) 



r, J/ 2 



ln(V27r) - 3/4 0.359933 



3/2 



Using the "robust" interpolation proposed by CioslowskP^ and the high- and low-density expansions ( 19 ) and ( 20 ) , 
the correlation energy can be approximated by 

3 

e c {r a )=t 2 Y,^{l-tf-\ (21) 

3=0 

with 

< - v/ "^" 1 . < 22 > 

2kr s 

and 

c = kr]o, ci =4fcr/ + /c 3/2 ?7i, (23) 

C2 = 5e + ei/A;, c 3 = ei, (24) 

where fc = 0.3083 is a scaling factor which is determined by a least-square fit of the DMC data given in Ref. |2H 

We disagree with the last comment made in Ref. |M1 which claims that this type of interpolation is not applicable to 
cases where the high- and low-density asymptotic expansions pertain to de facto different states, e.g. the 3D UEG. We 
claim that the non-applicability of such a interpolation is only due to the presence of logarithmic terms in the 2D and 
3D UEGs. However, in our case, the ID UEG does not involve any non-analytical term. Thus, the methodology of 
Ref. EHI is obviously applicable in the present case. 

The results using the new correlation functional ( 21 ) are compared to the DMC calculations of Lee and DrummoncP^ 
and to our recently proposed correlation functional based on a ring geometry^ and the "interaction-strength inter- 
polation" (ISI) formula proposed by Seidl and co-workers! 67 * 68 ! The results are gathered in Table |TT| and depicted in 
Fig. [2] As one can see, our new correlation functional reduces the error of the ISI functional by a factor two (expect 
for r s = 1 where the ISI functional is in better agreement with the DMC value). This overall improvement is expected 



because our new functional ( 21 ) contains one more information on the high-density limit than the correlation functional 
reported in Ref. [301 



IV. CONCLUSION 



In this Article, we have shown that the expression of the high-density correlation energy for the ID UEG is 
e c( r s) = —0.02742 + 0.00714r s + . . .. Combining these new results with the low-density correlation energy expansion 
£c(r s ) — — [ln(\/27r) — 3/4] rj 1 + 0.359933 r s 3 ^ 2 + . . . and the available diffusion Monte Carlo data, we have proposed a 
LDA correlation functional, which yields satisfactory estimates of the correlation energy at high, intermediate and low 
densities. 
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